The Abundance of KRAS and RAS Gene Mutations in Cancer.

Mutant forms of the RAS genes KRAS, NRAS, and HRAS are important and common drivers of cancer. Recently, two independent teams that integrated cancer genomics with cancer epidemiology estimated that approximately 15-20% of all human cancers harbor a mutation in one of these three RAS genes. These groups also estimate KRAS mutations occur in 11-14% of all human cancers. Although these estimates are lower than many commonly encountered values, these estimates continue to rank KRAS and the ensemble of RAS oncogenes among the most common genetic drivers of cancer across all forms of malignancy.

Decreasing alloimmunization-specific mortality in sickle cell disease in the United States: Cost-effectiveness of a shared transfusion resource.

Red blood cell alloimmunization and consequent delayed hemolytic transfusion reaction (DHTR) incidence and mortality in patients with sickle cell disease (SCD) are high. A shared transfusion resource has decreased both in other countries, while in the United States cost concerns persist. We conducted a Markov cohort simulation of a birth cohort of alloimmunized patients with SCD to estimate lifetime DHTR incidence, DHTR-specific mortality, quality-adjusted life expectancy (QALE), and costs with the implementation of a shared transfusion resource to identify antibody history versus without (i.e., status quo). We conducted our analysis using a lifetime analytic time horizon and from a United States health system perspective. Implementation of shared transfusion resource projects to decrease cumulative DHTR-specific mortality by 26% for alloimmunized patients with SCD in the United States, relative to the status quo. For an average patient population of 32 000, this intervention would generate a discounted increment of 4000 QALYs at an incremental discounted cost of $0.3 billion, resulting in an incremental cost-effectiveness ratio of $75 600/QALY [95% credible interval $70 200-81 400/QALY]. The results are most sensitive to the baseline lifetime medical expenditure of patients with SCD. Alloantibody data exchange is cost-effective in 100% of 10 000 Monte Carlo simulations. The resource would theoretically need a minimum patient population of 1819 patients or cost no more than $5.29 million annually to be cost-effective. By reducing DHTR-specific mortality, a shared transfusion resource in the United States projects to be a life-saving and cost-effective intervention for patients with SCD in the United States.

Theoretical analysis reveals a role for RAF conformational autoinhibition in paradoxical activation.

RAF kinase inhibitors can, under certain conditions, increase RAF kinase signaling. This process, which is commonly referred to as 'paradoxical activation' (PA), is incompletely understood. We use mathematical and computational modeling to investigate PA and derive rigorous analytical expressions that illuminate the underlying mechanism of this complex phenomenon. We find that conformational autoinhibition modulation by a RAF inhibitor could be sufficient to create PA. We find that experimental RAF inhibitor drug dose-response data that characterize PA across different types of RAF inhibitors are best explained by a model that includes RAF inhibitor modulation of three properties: conformational autoinhibition, dimer affinity, and drug binding within the dimer (i.e., negative cooperativity). Overall, this work establishes conformational autoinhibition as a robust mechanism for RAF inhibitor-driven PA based solely on equilibrium dynamics of canonical interactions that comprise RAF signaling and inhibition.

Computational Random Mutagenesis to Investigate RAS Mutant Signaling.

This chapter describes how mathematical models can be used to investigate the possible range of behaviors for mutant forms of a protein. A mathematical model of the RAS signaling network that has previously been developed and applied to specific RAS mutants will be adapted for the process of computational random mutagenesis. By using this model to computationally investigate the range of RAS signaling outputs that would be anticipated over a wide range of the relevant parameter space, one can gain intuition about the types of behaviors that would be demonstrated by biological RAS mutants.

Making inroads to precision medicine for the treatment of autoimmune diseases: Harnessing genomic studies to better diagnose and treat complex disorders.

Precision Medicine is an emerging approach for disease treatment and prevention that takes into account individual variability in genes, environment, and lifestyle. Autoimmune diseases are those in which the body's natural defense system loses discriminating power between its own cells and foreign cells, causing the body to mistakenly attack healthy tissues. These conditions are very heterogeneous in their presentation and therefore difficult to diagnose and treat. Achieving precision medicine in autoimmune diseases has been challenging due to the complex etiologies of these conditions, involving an interplay between genetic, epigenetic, and environmental factors. However, recent technological and computational advances in molecular profiling have helped identify patient subtypes and molecular pathways which can be used to improve diagnostics and therapeutics. This review discusses the current understanding of the disease mechanisms, heterogeneity, and pathogenic autoantigens in autoimmune diseases gained from genomic and transcriptomic studies and highlights how these findings can be applied to better understand disease heterogeneity in the context of disease diagnostics and therapeutics.

Co-targeting KRAS G12C and EGFR reduces both mutant and wild-type RAS-GTP.

The combination of KRAS G12C inhibitors with EGFR inhibitors has reproducibly been shown to be beneficial. Here, we identify another benefit of this combination: it effectively inhibits both wild-type and mutant RAS. We believe that targeting both mutant and wild-type RAS helps explain why this combination of inhibitors is effective.

How often is each gene mutated within the cancer patient population?

Genome sequenced samples from cancer patients helped identify roles of different mutation types and enabled targeted therapy development. However, critical questions like what are the gene mutation rates among the patients? or what genes are most commonly mutated, pan-cancer? have only been recently answered. Here, we highlight this recent advance.

Identification of RAS mutant biomarkers for EGFR inhibitor sensitivity using a systems biochemical approach.

Mutations can be important biomarkers that influence the selection of specific cancer treatments. We recently combined mathematical modeling of RAS signaling network biochemistry with experimental cancer cell biology to determine why KRAS G13D is a biomarker for sensitivity to epidermal growth factor receptor (EGFR)-targeted therapies. The critical mechanistic difference between KRAS G13D and the other most common KRAS mutants is impaired binding to tumor suppressor Neurofibromin (NF1). Here, we hypothesize that impaired binding to NF1 is a "biophysical biomarker" that defines other RAS mutations that retain therapeutic sensitivity to EGFR inhibition. Both computational and experimental investigations support our hypothesis. By screening RAS mutations for this biophysical characteristic, we identify 10 additional RAS mutations that appear to be biomarkers for sensitivity to EGFR inhibition. Altogether, this work suggests that personalized medicine may benefit from migrating from gene-based and allele-based biomarker strategies to biomarkers based on biophysically defined subsets of mutations.

Cancer gene mutation frequencies for the U.S. population.

Mutations play a fundamental role in the development of cancer, and many create targetable vulnerabilities. There are both public health and basic science benefits from the determination of the proportion of all cancer cases within a population that include a mutant form of a gene. Here, we provide the first such estimates by combining genomic and epidemiological data. We estimate KRAS is mutated in only 11% of all cancers, which is less than PIK3CA (13%) and marginally higher than BRAF (8%). TP53 is the most commonly mutated gene (35%), and KMT2C, KMT2D, and ARID1A are among the ten most commonly mutated driver genes, highlighting the role of epigenetic dysregulation in cancer. Analysis of major cancer subclassifications highlighted varying dependencies upon individual cancer drivers. Overall, we find that cancer genetics is less dominated by high-frequency, high-profile cancer driver genes than studies limited to a subset of cancer types have suggested.

Mathematical Modeling to Study KRAS Mutant-Specific Responses to Pathway Inhibition.

This chapter will describe how mathematical modeling allows the RAS pathway to be studied with computational experiments. The mathematical model utilized simulates the biochemical reactions that regulate RAS signaling. This type of model incorporates knowledge of reaction mechanisms, including measured quantitative parameters that characterize these reactions for both wild-type and mutant RAS proteins. For an illustrative example, this chapter focuses on how modeling provided new insights that helped solve a problem that challenged the RAS community for nearly a decade: why do colorectal cancers with the KRAS G13D mutation, but not the other common KRAS mutations, benefit from EGFR inhibition? The methods described include computational dose-response experiments and the use of "computational chimeric" RAS mutants.

KRAS-Mutated, Estrogen Receptor-Positive Low-Grade Serous Ovarian Cancer: Unraveling an Exceptional Response Mystery.

We report on a woman with aggressive estrogen receptor-positive, KRAS-mutated ovarian cancer who achieved a remarkable response to combination therapy with the MEK inhibitor (trametinib) and the aromatase inhibitor (letrozole), even though the disease had failed to respond to a combination of a PI3K inhibitor and different MEK inhibitor, as well as to trametinib and the estrogen modulator, tamoxifen, and to letrozole by itself. The mechanism of action for exceptional response was elucidated by in vitro experiments that demonstrated that the fact that tamoxifen can have an agonistic effect in addition to antagonist activity, whereas letrozole results only in estrogen depletion was crucial to the response achieved when letrozole was combined with an MEK inhibitor. Our current observations indicate that subtle variations in mechanisms of action of outwardly similar regimens may have a major impact on outcome and that such translational knowledge is critical for optimizing a precision medicine strategy. KEY POINTS: This report describes the remarkable response of a patient with KRAS-mutated, estrogen receptor-positive low-grade serous ovarian cancer treated with trametinib (MEK inhibitor) and letrozole (aromatase inhibitor), despite prior progression on similar agents including tamoxifen (estrogen modulator). In vitro investigation revealed that tamoxifen can have agonistic in addition to antagonistic effects, which could be the reason for the patient not responding to the combination of trametinib and tamoxifen. The current observations suggest that drugs with different mechanisms of action targeting the same receptor may have markedly different anticancer activity when used in combinations.

Discernment between candidate mechanisms for KRAS G13D colorectal cancer sensitivity to EGFR inhibitors.

Phase three clinical trial evidence suggests that colorectal cancers with the KRAS G13D mutation may benefit from EGFR inhibitors, like cetuximab, in contrast to the other most common KRAS mutations. A mechanism to explain why this mutation behaves differently from other KRAS mutations had long been lacking. Two recent studies have reproduced KRAS G13D specific sensitivity to cetuximab in cellular models, and both have implicated the tumor suppressor NF1 as a critical variable in determining sensitivity and resistance. One study proposes a mechanism that focuses on the inhibition of active, GTP-bound wild-type RAS, which is proposed to occur to a greater extent in KRAS G13D tumors due to the inability of KRAS G13D to bind NF1 well. The other study suggests NF1 can convert GTP-bound KRAS G13D to inactive, GDP-bound KRAS G13D. Here, we report an inability to reproduce cellular and biophysical studies that suggested NF1 has strong GTPase activity on KRAS G13D. We also report additional data that further suggests only WT RAS-GTP levels are reduced with EGFR inhibition and that KRAS G13D is impaired in binding to NF1. These new experiments further support a mechanism in which cetuximab inhibits wild-type (HRAS and NRAS) signals in KRAS G13D colorectal cancers. Video Abstract.

Real-world data from a molecular tumor board demonstrates improved outcomes with a precision N-of-One strategy.

Next-generation sequencing (NGS) can identify novel cancer targets. However, interpreting the molecular findings and accessing drugs/clinical trials is challenging. Furthermore, many tumors show resistance to monotherapies. To implement a precision strategy, we initiated a multidisciplinary (basic/translational/clinical investigators, bioinformaticians, geneticists, and physicians from multiple specialties) molecular tumor board (MTB), which included a project manager to facilitate obtaining clinical-grade biomarkers (blood/tissue NGS, specific immunohistochemistry/RNA expression including for immune-biomarkers, per physician discretion) and medication-acquisition specialists/clinical trial coordinators/navigators to assist with medication access. The MTB comprehensively reviewed patient characteristics to develop N-of-One treatments implemented by the treating physician's direction under the auspices of a master protocol. Overall, 265/429 therapy-evaluable patients (62%) were matched to ≥1 recommended drug. Eighty-six patients (20%) matched to all drugs recommended by MTB, including combinatorial approaches, while 38% received physician's choice regimen, generally with unmatched approach/low degree of matching. Our results show that patients who receive MTB-recommended regimens (versus physician choice) have significantly longer progression-free (PFS) and overall survival (OS), and are better matched to therapy. High (≥50%) versus low (<50%) Matching Score therapy (roughly reflecting therapy matched to ≥50% versus <50% of alterations) independently correlates with longer PFS (hazard ratio [HR], 0.63; 95% confidence interval [CI], 0.50-0.80; P < 0.001) and OS (HR, 0.67; 95% CI, 0.50-0.90; P = 0.007) and higher stable disease ≥6 months/partial/complete remission rate (52.1% versus 30.4% P < 0.001) (all multivariate). In conclusion, patients who receive MTB-based therapy are better matched to their genomic alterations, and the degree of matching is an independent predictor of improved oncologic outcomes including survival.

The Interpretation of SARS-CoV-2 Diagnostic Tests.

Testing for SARS-CoV-2 has attracted a tremendous amount of attention as a tool to manage the ongoing COVID-19 pandemic. Although diagnostic laboratory testing is used ubiquitously by physicians and encountered regularly by individuals receiving medical care, several aspects of test interpretation are incompletely understood by medical communities and the general population, creating a significant challenge in minimizing the damage caused by disease spread through informed decision making and proper testing utilization. Here, general principles of test interpretation are reviewed and applied to specific examples, such as whether asymptomatic individuals should be tested, what it means to test positive (or negative), and how to interpret tests for "immunity passports." Unexpectedly, the answers seem to run contrary to many of the popular narratives about testing as a tool for managing COVID-19. Although testing is an important and essential part of managing diseases such as COVID-19, improper utilization can have unintended negative consequences.

A mechanism for the response of KRASG13D expressing colorectal cancers to EGFR inhibitors.

Previous analysis of Phase 3 clinical trial data for colorectal cancer patients treated with cetuximab revealed that patients harboring a KRAS mutation did not benefit from treatment. This finding set the stage for one of the first examples of cancer personalized medicine. Confusingly, patients with a Glycine to Aspartic Acid mutation at amino acid 13 of KRAS (KRASG13D) appeared to respond positively to cetuximab, suggesting this mutation is an exception to the rule that KRAS mutations confer resistance to Epidermal Growth Factor Receptor (EGFR) inhibitors. Oncologists have stated that the mechanism that explains why the KRASG13D mutation is an exception should be identified before KRASG13D colorectal cancer patients should be treated differently. We have recently elucidated this mechanism using mathematical modeling of the KRAS biochemical system coupled with experimental biology. The mechanism we revealed involves a cetuximab-mediated reduction in HRAS and NRAS signaling within KRASG13D cancer cells, owing to impaired binding of KRASG13D to the tumor suppressor, Neurofibromin (NF1).

A systems mechanism for KRAS mutant allele-specific responses to targeted therapy.

Cancer treatment decisions are increasingly guided by which specific genes are mutated within each patient's tumor. For example, agents inhibiting the epidermal growth factor receptor (EGFR) benefit many colorectal cancer (CRC) patients, with the general exception of those whose tumor includes a KRAS mutation. However, among the various KRAS mutations, that which encodes the G13D mutant protein (KRASG13D) behaves differently; for unknown reasons, KRASG13D CRC patients benefit from the EGFR-blocking antibody cetuximab. Controversy surrounds this observation, because it contradicts the well-established mechanisms of EGFR signaling with regard to RAS mutations. Here, we identified a systems-level, mechanistic explanation for why KRASG13D cancers respond to EGFR inhibition. A computational model of RAS signaling revealed that the biophysical differences between the three most common KRAS mutants were sufficient to generate different sensitivities to EGFR inhibition. Integrated computation with experimentation then revealed a nonintuitive, mutant-specific dependency of wild-type RAS activation by EGFR that is determined by the interaction strength between KRAS and the tumor suppressor neurofibromin (NF1). KRAS mutants that strongly interacted with and competitively inhibited NF1 drove wild-type RAS activation in an EGFR-independent manner, whereas KRASG13D weakly interacted with and could not competitively inhibit NF1 and, thus, KRASG13D cells remained dependent on EGFR for wild-type RAS activity. Overall, our work demonstrates how systems approaches enable mechanism-based inference in genomic medicine and can help identify patients for selective therapeutic strategies.

Modeling cell line-specific recruitment of signaling proteins to the insulin-like growth factor 1 receptor.

Receptor tyrosine kinases (RTKs) typically contain multiple autophosphorylation sites in their cytoplasmic domains. Once activated, these autophosphorylation sites can recruit downstream signaling proteins containing Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains, which recognize phosphotyrosine-containing short linear motifs (SLiMs). These domains and SLiMs have polyspecific or promiscuous binding activities. Thus, multiple signaling proteins may compete for binding to a common SLiM and vice versa. To investigate the effects of competition on RTK signaling, we used a rule-based modeling approach to develop and analyze models for ligand-induced recruitment of SH2/PTB domain-containing proteins to autophosphorylation sites in the insulin-like growth factor 1 (IGF1) receptor (IGF1R). Models were parameterized using published datasets reporting protein copy numbers and site-specific binding affinities. Simulations were facilitated by a novel application of model restructuration, to reduce redundancy in rule-derived equations. We compare predictions obtained via numerical simulation of the model to those obtained through simple prediction methods, such as through an analytical approximation, or ranking by copy number and/or KD value, and find that the simple methods are unable to recapitulate the predictions of numerical simulations. We created 45 cell line-specific models that demonstrate how early events in IGF1R signaling depend on the protein abundance profile of a cell. Simulations, facilitated by model restructuration, identified pairs of IGF1R binding partners that are recruited in anti-correlated and correlated fashions, despite no inclusion of cooperativity in our models. This work shows that the outcome of competition depends on the physicochemical parameters that characterize pairwise interactions, as well as network properties, including network connectivity and the relative abundances of competitors.

Proceedings of the fifth international RASopathies symposium: When development and cancer intersect.

This report summarizes and highlights the fifth International RASopathies Symposium: When Development and Cancer Intersect, held in Orlando, Florida in July 2017. The RASopathies comprise a recognizable pattern of malformation syndromes that are caused by germ line mutations in genes that encode components of the RAS/mitogen-activated protein kinase (MAPK) pathway. Because of their common underlying pathogenetic etiology, there is significant overlap in their phenotypic features, which includes craniofacial dysmorphology, cardiac, cutaneous, musculoskeletal, gastrointestinal and ocular abnormalities, neurological and neurocognitive issues, and a predisposition to cancer. The RAS pathway is a well-known oncogenic pathway that is commonly found to be activated in somatic malignancies. As in somatic cancers, the RASopathies can be caused by various pathogenetic mechanisms that ultimately impact or alter the normal function and regulation of the MAPK pathway. As such, the RASopathies represent an excellent model of study to explore the intersection of the effects of dysregulation and its consequence in both development and oncogenesis.

Quantitative Systems Pharmacology Analysis of KRAS G12C Covalent Inhibitors.

KRAS has proven difficult to target pharmacologically. Two strategies have recently been described for covalently targeting the most common KRAS mutant in lung cancer, KRAS G12C. Previously, we developed a computational model of the processes that regulate Ras activation. Here, we use this model to investigate KRAS G12C covalent inhibitors. We updated the model to include Ras protein turnover, and validation demonstrates that our model performs well in areas of G12C targeting where conventional wisdom struggles. We then used the model to investigate possible strategies to improve KRAS G12C inhibitors and identified GEF loading as a mechanism that could improve efficacy. Our simulations also found resistance-promoting mutations may reverse which class of KRAS G12C inhibitor inhibits the system better, suggesting that there may be value to pursuing both types of KRAS G12C inhibitors. Overall, this work demonstrates areas in which systems biology approaches can inform Ras drug development.

A phase 2 trial of personalized cytotoxic therapy based on tumor immunohistochemistry in previously treated metastatic pancreatic cancer patients.

BACKGROUND: The choice of a regimen in metastatic pancreatic cancer patients following progression on 1st line therapy is empiric and outcomes are unsatisfactory. This phase II study was performed to evaluate the efficacy of therapy selected by immunohistochemistry (IHC) in these patients following progression after one or more therapies. METHODS: Eligible patients underwent a percutaneous biopsy of a metastatic lesion and treatment selection was determined by IHC. The study required 35 evaluable patients (power of 86%) for detecting a true 1-year survival rate of >20%. RESULTS: A tumor biopsy was performed in 48 of 49 accrued patients. Study therapy was not given (n=13) either due to insufficient tumor on biopsy (n=8) or due to worsening cancer related symptoms after biopsy (n=5). The demographics of evaluable patients (n=35) are male/female (59%/41%), with age range 34-78 years (median 63 years). Patients had 1-6 prior regimens (median of 2). The most common IHC targets were topoisomerase 1 or 2, thymidylate synthase, excision repair cross-complementation group 1 protein (ERCC1), and osteonectin secreted protein acidic and rich in cysteine (SPARC). Commercially available treatment regimens prescribed included FOLFIRI, FOLFOX, irinotecan, and doxorubicin. The response (RECIST) was 9%, the median survival was 5.6 months (94% CI, 3.8-8.2), and the 1-year survival was 20% (95% CI, 7-33%). CONCLUSIONS: In all patients, IHC assays resulted in identification of at least two targets for therapy and a non-cross resistant regimen could be prescribed for therapy with evidence of some benefit. An IHC based treatment strategy is feasible and needs validation in larger studies.